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Abstract 

Omnivory is defined as feeding on more than one trophic level. An 
example of this is the so-called intraguild predation (IG) which includes a 
predator and its prey that share a common resource. IG predation mod¬ 
els are known to exhibit interesting dynamics including chaos. This work 
considers a three-species food web model with omnivory, where the inter¬ 
actions between the basal resource, the IG prey, and the IG predator are 
of Lotka-Volterra type. In the absence of predation, the basal resource 
follows a delayed logistic equation or popularly known as Hutchinson’s 
equation. Gonditions for the existence, stability, and bifurcations of all 
non-negative equilibrium solutions are given using the delay time as pa¬ 
rameter. Results are illustrated using numerical bifurcation analysis. 
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1 Introduction 

Omnivory, as defined in occurs when a population feeds on resources at 
more than one trophic level. For example, a species feeding on its prey’s resource 
is called omnivorous. This particular type of tri-trophic community module is 
called intraguild (IG) predation |I6j , and was shown to be quite common in na¬ 
ture [1]. This three-species IG predation model includes top and intermediate 
predators termed as IG predator and IG prey, respectively, and a basal resource. 
The IG predator depends completely both on the IG prey and the basal resource 
for its sustenance, while the IG prey depends solely on the basal resource. A 
study of a model of IG predation with non-linear functional responses in m 
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showed that omnivory stabilizes and enhances persistence of the three-species 
food web. However, a model of Lotka-Volterra type with linear functional re¬ 
sponses considered in [S] showed that IG predation could have a destabilizing 
effect, and a criterion for co-existence of all three species is that the IG prey 
must be superior than the IG predator in competing for the shared basal re¬ 
source. Lotka-Volterra IG predation models are known to exhibit interesting 
dynamics such as limit cycles [ 8 ], bistabity [ 10 ], and chaos nain]. 

In this paper, we consider a Lotka-Volterra IG predation model where, in the 
absence of the IG predator and the IG prey, the basal resource grows according 
to the delayed logistic equation or more commonly known as the Hutchinson’s 
equation l^. It should be noted that there are several other ways of formulat¬ 
ing the delayed logistic equation such as in | 2 ], and in this paper, we use the 
so-called classical delayed logistic equation. Our delayed model generalizes the 
Lotka-Volterra IG predation models of [Sj and m which are basically systems 
of ODEs. 

This paper is organized as follows. In Section 2, we introduce the delayed 
Lotka-Volterra IG predation model and discuss the existence of its equilibrium 
solutions. In Section 3, we give the main results of this paper. These are 
the conditions for stability and bifurcations of all non-negative equilibria using 
the delay time as parameter. In Section 4, we use numerical continuation and 
bifurcation analysis to illustrate our results on the effects of the delay time to 
our Lotka-Volterra IG predation model. We then end by giving a summary and 
the thoughts of this paper. 

2 The Model 

We consider the following delayed Lotka-Volterra IG predation model 

x{t) = [ao - aix{t - t) - a2y{t) - a3z{t)]x(t), 

y{t) = [-bo + bix{t) - b3z{t)]y{t), ( 1 ) 

z{t) = [-co +cix{t) + C2y{t)]z{t), 

where x{t), y{t), and z{t) are the densities at time t of the basal resource, IG 
prey, and IG predator, respectively. In the absence of the IG predator and 
the IG prey, the basal resource grows according to the delay logistic equation 
|3|[T8]. When r = 0, we recover the ODEs model in [19]. We refer to [19] for 
the description of the rest of the parameters in ([1]). Moreover, we use an initial 
condition (x[t),y(t), z(t)) = {xo,yo,zo) for t € [— t , 0 ] and where xo,yo,zo > 0 . 

An important characteristic of the IG predation is that is it a mixture of com¬ 
munity modules such as competition and predation [8j|9j[T6]. When 03 = ci = 0, 
we obtain a food chain while if 63 = C 2 = 0 , we obtain exploitative competition 
(or shared resources) where two predators, in our case the IG predator and the 
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IG prey, share a common resource and the IG predator does not feed on the IG 
prey. Meanwhile, if 02 = = 0 , we obtain apparent competition (or shared 

predation), which in our model, the IG predator feeds on both the IG prey and 
the basal resource but there is no predation on the basal resource by the IG 
prey. In this case, the IG prey will go extinct. 

2.1 Equilibrium Solutions 

System © has five possible non-negative equilibrium solutions: Eq = (0,0,0), 
El = (K, 0 , 0 ) where K = ao/oi, E2 = (A,B, 0 ) where A = 60/61 and B = 
(ao6i — 0160)70261, E3 = (C, 0 , D) where C = co/ci and D = (ooCi — aiCo)/a3Ci, 
and the positive equilibrium solution E/i = {P/S,Q/S, R/S) where 


P = 006302 - 0263C0 -I- 036002, 

Q = -0063O1 -I- 016300 - 036001 -I- 036100, 

R = 0061O2 — 0160O2-I-0260O1 — 026100, 

S = 016302 — 026301 -l- 036102. 


Since all parameters of system 0 are positive, we have K > 0, A > 0, and 
C > 0. Thus, Eq and Ei always exist. For E 2 and Eq to exist, we require B > 0 
and D > 0, respectively, or equivalently A < K and C < K, respectively. Now, 
if S is positive (resp. negative), then each of P, Q, and R must also be positive 
(resp. negative) for E 4 to be a positive equilibrium. The following theorem 
summarizes these results. 

Theorem 2.1 For system the equilibrium solutions Eq = (0,0,0) and 

El = (A', 0, 0) always exist, while E 2 and Eq exist provided agbi > 0160 and 
oooi > oioo, respectively, or equivalently, A < K and C < K, respectively. The 
positive equilibrium solution E 4 exists if P/S, Q/S, and R/S are all positive. 

2.2 Local Stability of Equilibria 

Let X{t) = [x{t),y{t), z{t)]'^. Then, the linear system corresponding to ([T]) 
around an equilibrium solution A* = (a:*,!/*, 2 *) isX{t) = MQX(t) + MiX{t—T) 
where 


oq — aicc* — a2y* — oqz^, 


—a 2 X^, 

bo + bix^. — 632 ;* 
C2Z* 


— OqX^ 


Mq 


biy* 

ClZ* 


Co -I- cia:* -I- C2y. 


-hy* 


and 


Ml 


—aiXt, 0 0 

0 0 0 

0 0 0 


and with corresponding characteristic equation 


det(A/ - Mo - Mie"^^) = 0 . 


( 2 ) 
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If all roots of ([2]) have negative real part, then the equilibrium solution is 
locally asymptotically stable [18] . We can think of the roots of (l2|) as continuous 
functions of the delay time r, that is, A = A(r). When r = 0, ([2]) is a polynomial 
equation. In this case, the well-known Routh-Hurwitz criterion can then be uti¬ 
lize to provide stability conditions for E^:. As r is increased from zero, Corollary 
2.4 of [T7] tells us that the sum of the orders of the roots of (I2|) in the open right 
half-plane can change only if a zero appears on or crosses the imaginary axis. 
That is, stability switch occurs at a critical value t = tq where A(to) is either 
zero or purely imaginary. The transversality condition d{Re X) /dT\.r=To > 0 
implies that the eigenvalues cross from left to right. Hence, if A* is stable at 
r = 0, then, as r is increased, it loses its stabilty at r = tq. Thus, A* is locally 
asymptotically stable when r S (0 ,to). If there are no roots of (l2|) that cross 
the imaginary axis, then there are no stability switches and in this case, we have 
absolute stability |3]. That is, A* remains stable for all delay time r > 0. 

3 Main Results 

In the following, we give conditions for the stability of each equilibrium solution 
of system (HD when T > 0 using the technique mentioned above. When r = 0, 
the stability analysis of each of the equilibrium solution of system (HD can be 
found in |141 Appendix A] and in [^ but in slightly different form. We mention 
them here for completeness. 

At Eg, the characteristic equation ([2|) becomes (A — ao)(A -I- bo)(A + cg) = 0 
whose roots are ag >0, —bg < and — cq < 0. Thus, we have the following result. 

Theorem 3.1 The equilibrium solution Eg = (0,0,0) of system (Qp is a local 
saddle point and is unstable for all delay time r > 0. 

At Ai, the characteristic equation (l 2 |) becomes 

(A — bxa-iBjai) (A — ciagDjai) (A -I- age~^'^) = 0. (3) 

When r = 0, dH) has roots —ag < 0, &ia 2 H/ai, and cxagDjax. In this case, 
we see that Ei is locally asymptotically stable provided both B and D are 
negative, or equivalently A > K and C > K. Recall from Theorem 12.11 that 
these conditions imply that both E2 and Eg do not exist. Suppose now that 
r > 0, and A > K and C > K. It is known that all roots of A -I- = 0 

have negative real part if and only if agr < 7 r /2 (see for example [T2J pp.70]). 
Thus, El is locally asymptotically stable if and only if r < 7 r/ 2 ao. 

Theorem 3.2 Suppose that m satisfies aibg > agbi and aicg > agCi, or 
equivalently A > K and C > K, respectively. Then, the equilibrium solution 
El = (if, 0,0) is locally asymptotically stable if and only if t € (0 ,Tc) where 
Tc = 7 r/ 2 ao. If T = Tc, then ([I]) undergoes a Hopf bifurcation at Ei. 
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Example 1. The set of parameter values 


(flo, Oi, 02, 03, 6oi bi, 63, Co, Cl, C2) 


(1,0.5,1,0.6,0.75,0.25,0.5,0.5,0.15,0.3) 


satisfies the assumptions in Theorem 13.21 and gives Tc = 7 r /2 (approx. 1.5708). 
Figure [T] illustrates the stability switch at r = 7 r /2 using the initial condition 
(x(t), y{t),z{t)) = ( 2 , 1 , 1 ) for t G [—r, 0 ]. Ei = ( 2 , 0 , 0 ) is stable for r G ( 0 , Tc). 


X = 1.50 




time t 


x = 1.60 




0 50 100 150 



0 50 100 150 


time t 


Figure 1: Stability switch occurs at a Hopf bifurcation point where r = tt/2. 
El is stable when t <7112 (Left) and is unstable when r > 7 r /2 (Right). 


At £' 2 , the characteristic equation ([ 2 |) becomes 

(A-R/a2&i)(A2+6Ae-^^ + c) =0 (4) 

where b = aiA and c = a 2 boB. Since the existence of £2 requires both A and 
B to be positive, we have 6 > 0 and c > 0. First, consider the case when r = 0. 
Since b > 0 and c > 0, the Routh-Hurwitz criterion tells us that all roots of 
A^ + 6 A + c = 0 have negative real parts. Thus, £2 is locally asymptotically 
stable if and only if £ < 0. Suppose now that r > 0 and £ < 0. Since both c and 
£ are non-zero, then A = 0 is not a root of (|H). Meanwhile, notice that A = i/i, 
with y > 0, is a root of dH) if and only if y satisfies + ibye~'‘^'^ -b c = 0. 
Splitting into real and imaginary parts, we obtain 

hysin{yT) = y? — c and 6 /r cos(/iT) = 0. (5) 


Squaring each sides of equations in ([5|) and then adding them together, we get 


6 ^ -b 2 c ± \J (52 -p 2 c) ^ — 4 c^ 


( 6 ^ -b 2 c) /r^-bc^ = 0 , or equivalently, y? = ^ 

Since P -b 2 c > 0 and the discriminant (P -b 2 c)^ — 4c^ = IP (IP + 4c) > 0, we 
obtain two positive valnes of y, given by 


= I ^ {b'^ + 2 c) ± \J (62-b2c)^-4c2 I 


1/2 


( 6 ) 


5 





























with 0 < fi- < M+- Using we get cot fiT = 0. So that //r = {2k + l)7r/2 for 
fc = 0,1, 2,.... Now, corresponding to a given fi, define the increasing sequence 
Tk = {2k + l)7r/2/r for fc = 0,1, 2,.... The proof of the following lemma can be 
found in [TH pp. 74-75]. 

Lemma 3.3 Let A(t) be the root of satisfying A(r/c) = ifr. Then, 

= sign {2g,^ — 1? — 2c). 

T = Tk 


sign 


d{Re A) 
dr 


For ^ = yi±, define the sequence = {2k + l) 7 r/ 2 ^± for k = 0,1, 2,... cor¬ 
respondingly, and let Since 0 < /i_ < /i_|_, we have < r^_, and 

hence is the smallest amongst all r^. From ([51), we have 2^^ — 6 ^ — 2c > 0. 
Thus, at T = (or equivalently at ^ = /r+), the quantity is positive 

by Lemma 15751 This, together with the Hopf bifurcation theorem [T5], we have 
the following result. 

Theorem 3.4 Suppose that in system m, aobiC 2 + a 2 boCi < aibQC 2 + a 2 biCo 
or equivalently R < Q. Then, the equilibrium solution E 2 = {A,B,0) of (ifl is 
loeally asymptotieally stable whenever t S (0,t^^). If t = then system (Qp 
undergoes a Hopf bifurcation at E 2 ■ 


Example 2. Using the same set of parameters and initial condition as in 
Example 1 with bi changed to 0.50, the assumptions in Theorem l3.4l are satisfied. 
Figure [2| illustrates the stability switch at r = = 1.6573 approximately. 


X = 1.50 



0 50 100 150 



50 100 ISO 


time t 


x= 1.75 



0 50 100 150 


time t 


Figure 2: Stability switch occurs at a Hopf bifurcation where r = = 1.6573 

approximately. The equlibrium solution E 2 = (1.50,0.25,0.00) is stable when 
T < (Left) and is unstable when t > (Right). 

A similar condition for the stability of E 3 can be obtained using the same 
analysis used for E 2 , and is given in the following theorem. 
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Theorem 3.5 Suppose that in system m, aib^co + asbico < aobsci + a^boci 
or equivalently Q < 0. Let b = aiC, c = o^cqD, and 


(i^ + 2 c^ y 

Then, the equilibrium solution = (C, 0, D) of (QP is locally asymptotically 
stable whenever r G (0,Tjy^) where = tt/2vj^. If t = then system (Op 
undergoes a Hopf bifurcation at . 

Stability of the Positive Equilibrium 

At E 4 , the characteristic equation @ becomes 

X^ + aX + b+ (cA^ + d)e-^^ = 0 (7) 

where a = {PQa 2 bi + PRa^ci + QRb^C 2 )/S'^, b = PQR{a^biC 2 — a 2 b 3 Ci)/S^, 
c = aiP/S, and d = PQRaib 3 C 2 /. Note here that b + d = PQR/S^. Order- 
three quasi-polynomials with single delay, different to 0 . have been studied in 
[3 HI]. We follow a similar method used in [3 to analyze the distribution of the 
roots of 0 on the complex plane. When r = 0, then 0 reduces to 

-I- cA^ -I- oA -I- (6 -I- d) =0. ( 8 ) 

The Routh-Hurwitz criterion tells us that all roots of (jHj) have negative real part 
if and only if the following inequalities hold: a,c,b + d > 0 and ac — (b + d) > 0 . 
If 5 < 0, then P, Q, and R must all be negative for E 4 to be a positive equilib¬ 
rium. This implies that a > 0, c < 0, and 5 -I- d < 0, and thus, E 4 is unstable. 
If S' > 0, then each of P, Q, and R must be positive for E 4 to be a positive 
equilibrium. In this case, a, c, and b d are all positive. Now, notice that if 
in addition to the assumption that S > 0, we also have b < 0 (or equivalently 
026301 > 036102 ), then 00 — (6 -I- d) = aiP^{Qa 2 bi + Ra 3 Ci)/S^ — 6 > 0. Hence, 
for the case when t = 0, A 4 is locally asymptotically stable whenever S > 0 
and 6 < 0 . 

Suppose now that r > 0, and assume that S > 0 and 6 < 0. If A = 0 is a 
root of 0 , then 6 -I- d = 0. However, 6 -I- d 7 ^ 0 since 6 -|- d = PQR/S^ and none 
of P, Q, and R is equal to zero. Thus, A = 0 is not a root of (HI). If A = iio, 
with oj > 0, is a root d?]), then —iuj^ + iauj -|- 6 -|- {—cuP' P d)e“*“'^ = 0. Splitting 
into real and imaginary parts, we get 

(chP — d) cos LOT = 6 , and (ow^ — d) sinwr = — ow. (9) 

Squaring each sides of these equations and then adding corresponding sides gives 

w® -I- acj'^ + flop - 1-7 = 0 (10) 
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where a = —2a — (3 = a? 2cd, and 7 = 6 ^ — We claim that a < 0, 

(3 > 0, and 7 < 0. To see this, observe that since S' > 0, we have P,Q, R > 0 
and consequently, a, c, and d are all positive. Immediately, we see that a < 0 
and f3 > 0. Meanwhile, the assumptions S > 0 and 6 < 0 gives b + d > 0 and 
6 — d < 0, respectively. As a result, 7 = (d + d )(6 — d) < 0. We now show that 
m has at least one positive root. Let h{u) = + Pu + "f, and note that 

if the equation h{u) = 0 has a positive root u = uq, then pOI) has a positive 
root Wo = v^. Now, notice that h'{u) = 3u^ + 2au + /3 = 0 has two positive 
roots u± = {—a ± —‘3(3 )/3 with 0 < < u+. Since h{Q) = 7 < 0, 

d'(0) = /3 > 0, and h"(0) = 2a < 0, we know that the graph of h passes through 
the point (0,7) below the horizontal axis, and then increases in a concave down 
manner. The continuity of h and the fact that h increases without bound as 
u —>■ +00 guarantee that h has at least one positive root mq. That is, OT has a 
positive root wq = i/uo, and therefore ([71) has purely imaginary roots A = ±zwo. 
Noting that a local maximun of h occurs at it = it_ and a local minimun of h 
occurs at It = it+, we see that the graph of the cubic polynomial h is increasing 
on (0, it_) and on (it+, +00). Moreover, we either have 0 < itg < ii_ or mq > u+. 


If the equation h{u) = 0 has more than one positive roots, then it must 
have exactly three positive roots so that w is a simple root of (nni). Using the 
hrst equation in ([5]), for a given w, define its corresponding increasing sequence 

^ ( cuj^-d ) ■*" A: = 0,1, 2,.... Specifically, for a given w, we 


Tk = — 


cos 


get a corresponding tq. Among the three positive roots of the equation h{u) = 0, 
we then choose mq so that ui = ujq = ^/uq has corresponding tq that is smallest. 
This guarantees that a pair of purely imaginary eigenvalues will first occur at 

T = To. 


Lemma 3.6 Let A(r) be the root of & satisfying X{Tk) = iwq. Then, 


d{Re A) 


dr 


r=Tk 


> 0 . 


Proof. From ([7]), 


[(3A^ + a) + (2cA — (cA^ + d)r)e - -A(cA^ + d)e = 0. 

Consequently, 

/dA\ ^ 2cA + (3A^ + a)e^''’ r 
\dT) A(cA^ + d) A 

2c 3A^ + a T 

cA^ + d A(A^ + aX + b) X 
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since e'^’’/(cA^ + d) = —1/(A^ + aA + 6) using ([7]). Hence, 


sign 


f d {Re A) 'I 

I dr J 


X=iuJo 


sign -i 

^i?e| 

signj 

|i?e ■ 

signj 

■ 

signj 

r - 

LcWq 

sign <1 

f - 

[cWq 

sign j 

f - 

[cwl 


X—iojQ 


—3A^ — a 


X—iijQ 


—cuJq + d 


3wq — a 

Wq — awg + ibujo 


- d 


- d 


(3a;g - a){uj^ -auJo) 

(w^ - + (6a;o)2 

(3a;g -a)(a;g -g) ] 

- awo)^ + J 
(3^0 -q)(^o -q) 1 

(cc^g-d)2 / 


since (wq — awo)^ + 6 ^ = (cWq — from ([5]). Thus, 


sign 


d {Re A) 
dr 


^ X—iijjQ 


sign {—2 c(cwq — d) + (3a;o — g)(wo ~ ®)} 

sign {3a;o + (—4a — 2c^)a;o + (a^ + 2cd)} 
sign {3a;o + 2 q;Wq + /?} . 


Since h is increasing on the intervals (0,w_) and (w+,oo), and uq belongs to 
(0,M_) or to (u+,oo), we have h'{uo) = 3uq + 2aMo + ^5 > 0. This implies that 
3 wq + 2aujQ + /3 > 0 since wq = ^/uo■ Thus, 

= sign {3wo + 2 aa;o + /?} = +1 

X—iuJo 

and this completes the proof. 


sign 


d {Re A) 
dr 


By Lemma 13.61 


d{Re A) 
dr 


T — To 


> 0. This, together with the Hopf bifurcation 


theorem m give the following result. 


Theorem 3.7 Suppose that in system m, aib 3 C 2 + a3&iC2 > a2l'3Ci and 
a 2 & 3 Ci > a 3 biC 2 (or equivalently S > 0 and b < 0, respectively). Then, the 
positive equilibrium solution E 4 of w is locally asymptotically stable whenever 
T G (0,ro). If T = To, then system (1) undergoes a Hopf bifurcation at E 4 . 


Example 3. Using the same set of parameters as in Example 1 with 61 = 1 
and Cl = 0.42, the assumptions in Theorem 13.71 are satisfied. Using the initial 
condition {x{t),y{t), z{t)) = (0.78,0.58,0.06) for t G [—t, 0], Figure [3] illustrates 
the stability switch at r = tq = 1.7438 approximately. 
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Figure 3: Stability switch occurs at r = tq = 1.7438 (approx.). The positive 
equilibrium solution = (0.7778,0.5778,0.0556) is stable for r < tq (Left) and 
is unstable when t > tq (Right). 


4 Numerical Continuation 

Recall that for t = 0, E 4 is locally asymptotically stable whenever S > 0 and 
b < 0. A branch of equilibrium solutions can be obtained by following or con¬ 
tinuing this equilibrium solution in DDE-Biftool [6] using the delay time r as 
parameter. DDE-Biftool is a numerical continuation and bifurcation analysis 
tool developed by Engelborghs et al [7] . We use the same set of parameters and 
initial condition used in Example 3. In Figured the top panel shows the branch 
of equilibrium solutions (horizontal line) obtained in DDE-Biftool, where green 
and magenta represents the stable and unstable parts of the branch, respec¬ 
tively. A change of stability occurs at a Hopf bifurcation point marked with (*) 
where r = tq. Again, we use DDE-biftool to continue this Hopf point into a 
branch of periodic solutions. We obtained a stable branch of periodic solutions 
and this shown as the green curve in the top panel of Eigure|4l Here, the verti¬ 
cal axis gives a measure of the maximum value of the oscillation of the periodic 
solutions in the x(t) component. The middle and bottom panels corresponds 
to the same bifurcation diagram as the top panel but for the y(t) and the z(t) 
components. 

EigureS] illustrates the results of Theorem l3.71 that is, the stability switch at 
T = Tq and the occurrence of Hopf bifurcation. Numerical continuations showed 
that, as the positive equilibrium E^ loses its stability, the branch of periodic 
solutions that emerges from the Hopf point is stable. Biologically, this means 
that the three species still persist and the only difference now is that the three 
populations are oscillating. 
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Figure 4: Stable branch of periodic solutions emerges from the Hopf bifurcation 
point marked with (*) where r = tq. 


5 Conclusion 

We studied a delayed Lotka-Volterra IG predation model where, in the absence 
of predation, the basal resource grows according to the delayed logistic equa¬ 
tion. By first considering the case where r = 0, stability conditions for the all 
non-negative equilibria are established using the well-known Routh-Hurwitz cri¬ 
terion. For T > 0, we showed that stability switch occurs at a Hopf bifurcation 
where the stable equilibrium becomes unstable. Moreover, numerical continu¬ 
ation shows that a stable branch of periodic solutions emerges from the Hopf 
point as the positive equilibrium solution becomes unstable. This means that 
by increasing the delay time the positive equilibrium could become unstable. 
However, since the periodic orbit that will emerge is stable, all three species 
will still persist. 
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